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Abstract 

Using the methods of the "form factor program" exact expressions of all matrix 
elements are obtained for several operators of the quantum sine-Gordon model alias 
the massive Thirring model. A general formula is presented which provides form 
factors in terms of an integral representation. In particular charge-less operators as 
for example the current of the topological charge, the energy momentum tensor and 
all higher currents are considered. In the breather sector it is found the quantum 
sine-Gordon field equation holds with an exact relation between the "bare" mass 
and the normalized mass. Also a relation for the trace of the energy momentum is 
obtained. All results are compared with Feynman graph expansion and full agree- 
ment is found. 
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1 Introduction 

This work continues a previous investigation [[j] on exact form factors for the sine- Gordon 
alias the massive Thirring model. Coleman [TJJ had shown that these two models are 
equivalent on the quantum level. The corresponding classical models are defined by their 
Lagrangian's 

C SG = i^^+^(cos^-l) 
C MT = M)ip — \g j% , (f = VW) • 
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We do not use these classical Lagrangians and any quantization procedure to construct 
the quantum models. We have contact with the classical models only, when at the end 
we compare our exact results with Feynman graph expansions which are based on the 
Lagrangians. The 'form factor program' is part of the 'Bootstrap program for integrable 
quantum field theories in 1+1-dimensions'. This program classifies integrable quantum 
field theoretic models and in addition it provides their explicit exact solutions in term of 
all Wightman functions. These results are obtained in three steps: 



1. The S-matrix is calculated by means of general properties such as unitarity and 
crossing, the Yang-Baxter equations (which are a consequence of integrability) and 
the additional assumption of 'maximal analyticity'. This means that the two-particle 
S-matrix is an analytic function in the physical plane (of the Mandelstam variable 
(pi +P2) 2 ) and possesses only those poles there which are of physical origin. 

2. Generalized form factors which are matrix elements of local operators 

out (p' m ,...,p' 1 \0(x)\p 1 ,...,p n ) m 

are calculated by means of the S-matrix. More precisely, the equations (i) — (v) given 
below on page | are used as an input. These equations follow from LSZ-assumptions 
and again the additional assumption of 'maximal analyticity' (see also [Q]). 

3. The Wightman functions are obtained by inserting a complete set of intermediate 
states. In particular the two point function for a hermitian operator 0(x) reads 
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The on-shell program i.e. the exact determination of the scattering matrix was formu- 
lated in ||. Off-shell considerations were carried out in §1 and in ||, where the concept 
of a generalized form factor was introduced and consistency equations were formulated 
which are expected to be satisfied by these objects. Thereafter this approach was devel- 
oped further and studied in the context of several explicit models (see e.g. @[|). More 
recent papers on solitonic matrix elements in the sine-Gordon model are J7|, §]. There is a 
nice application || |Kj] of form factors in condensed matter physics. The one dimensional 
Mott insulators can be described in terms of the quantum sine-Gordon model. 

In the previous paper M an integral representation for general matrix elements of the 



fundamental fermi-field of the massive Thirring model has been proposed. In [IT], [12], 
we generalize this formula and investigate in particular charge-less local operators. The 
strategy is as follows: 

For a state of n particles of kind ol{ with momenta pi = m sinh 8i and a local operator 
0(x) the generalized form factor is defined by 

( I O(x) I ai(pi), . . . , a n ( Pn ) P = e- ix ^ + - + ^Oa(9) (1) 



1 For more references sec Q]. 
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for 6 1 > ■ ■ ■ > 9 n . The short notation a = (ax, . . . , a n ) and 9 = (9\, . . . , 6 n ) has been 
used. We make the Ansatz 

Oa(0)= [ d Zl --- f dz m h(9,z)p°(9,z)^^9,z) 

JC 6 _ JCg_ 

with the Bethe state ^a(9, z) defined by eq. ((5|) and the integration contours Cg of figure [3]. 
The scalar function h(9_,z_) is uniquely determined by the S-matrix whereas the scalar 'p- 
function' p°(9,z) depends on the operator. By means of the Ansatz we transform the 
properties (i) — (t>) of the co-vector valued function Og_(9) (see page |[) to properties 
(i') — (V) of the scalar function p°(9,z) which are easily solved. For example we obtain 
the p-functions for the local operatoif] Af \j/}ip\ (x) as 




In section |] we propose in addition the p-functions for J\f ^7 5 ^] (x), the current 
the energy momentum tensor T^ v {x) and the infinitely many higher conserved currents 
J^{x). The identification with the operators is made by comparing the exact results with 
Feynman graph expansions. Properties as charge, behavior under Lorentz transformations 
etc. will also become obvious. 



2 Recall of formulae 

In this section we recall some formulae which we shall need in the following sections to 
present our results. All this material can be found in [0 including the original references. 

2.1 The S-matrix 

The sine-Gordon model alias massive Thirring model describes the interaction of several 
types of particles: solitons, anti-solitons alias fermions and anti-fermions and a finite 
number of charge-less breathers, which may be considered as bound states of solitons and 
anti-solitons. Integrability of the model implies that the n-particle S-matrix factorizes 
into two particle S-matrices. In particular scattering conserves the number of particles 
and even their momenta. The two particle soliton-soliton amplitude a(9), the soliton 
anti-solitons forward and backward amplitudes b{9) and c(9) are given by p4f |3| 

ua\ sinhfl/z/ sinhivr/z/ 

¥ = ^7 m^ a #' c {9) = — — -r 7-o(0), 

smh(27r — 9)jv smh(27r — 9) jv 

f°° dt sinh|(l - v )t e 

a{9) = exp / -V sinht— . (2) 

J t sinh ±vt cosh ±t m 



2 The symbol M refers to normal products of local quantum fields. 
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The parameter 9 is the absolute value of the rapidity difference 9 — \ 9\ — 92 | where 9i are 
the rapidities of the particles given by the momenta Pi = Msinh#j. The parameter v is 
related to the sine-Gordon and massive Thirring model coupling constant by 

(3 2 

v 



where the second equality is due to Coleman M. 

We list some general properties of the two-particle S-matrix. As usual in this context 
we use in the notation 

v l...n e yl...n ^V X ®---®V n 

for a vector in a tensor product space. The vector components are denoted by v ctl '" (Xn . A 
linear operator connecting two such spaces with matrix elements Aa\..2n is denoted by 

A i:.n ' ■ yl '" n ^ 1 '"' n ' 

where we omit the upper indices if they are obvious. All vector spaces Vi are isomorphic 
to a space V whose basis vectors are label all kinds of particles (here solitons and anti- 
solitons, i.e. V = C 2 ). An S-matrix as acts nontrivial only on the factors V, ® Vj. 

The physical S-matrix in the formulas above is given for positive values of the rapidity 
parameter 9. For later convenience we will also consider an auxiliary matrix S(9i,9 2 ) 
regarded as a function depending on the individual rapidities of both particles 9 1 ,9 2 or 
some times also on the difference 9\ — 9 2 

6 in n\_o (a a \ _ / { (J S) 12 {\9 1 - 9 2 \) for 9 1 > 9 2 

b 12 {9 1 ,9 2 )-b 12 {9 l -9 2 )-^ (Sa) 2 l{\9 1 -9 2 \) for 9 1 < 9 2 

with a taking into account the statistics of the particles. It is a diagonal matrix o\ 2 with 
entries —1 if both particles are fermions and +1 otherwise. The matrix S(9x,9 2 ) is an 
analytic function in terms of both variables 9\ and 9 2 . The auxiliary matrix S12 acts 
nontrivial on the factors V\ <8> V 2 and in addition exchanges these factors, e.g. 



S 12 (9) : V 1 ® V 2 -> V 2 ® V x 



It may be depicted as 



S\ 2 (9i, 9 2 




Here and in the following we associate a rapidity variable ^ 6 C to each space Vi which 
is graphically represented by a line labeled by 9i or simply by i. In terms of the auxiliary 
S-matrix the Yang-Baxter equation has the general form 

5*12(^12) Si 3 (9i 3 ) S 2S (9 23 ) = S 23 (9 23 ) S ri (9i 3 ) Si 2 (9i 2 ) 
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which graphically simply reads 





Unitarity and crossing may be written and depicted as 



S , 2 l(^2l)'S'l2(^12) = 1 '■ 



12 12 



(3) 



£12(01 - 62) = C 22 5 2 i (02 + wr - 61) C22 = C 11 S 2l (9 2 - {B t - in) C 11 (4) 






1 2 



1 2 



where C 11 and Cji are charge conjugation matrices. For the sine-Gordon model the 
matrix elements are C a/3 = C a p = 5 af 3 where (3 denotes the anti-particle of (3. We have 
introduced the graphical rule, that a line changing the "time direction" also interchanges 
particles and anti-particles and changes the rapidity as 9 — > 9 ± in, as follows 



9 — in 



iaf3 



a 13 
9\ \9 + m 



a (3 



Similar crossing relations will be used below to investigate the properties of form factors. 
Finally we denote a property of the two-particle S-matrix 



which turns out to be true for all examples. This means that S for zero momentum 
difference is equal to minus the permutation operator. 



2.2 Form factors 

For a state of n particles of kind aij with momenta Pi and a local operator 0(x) we 
define the form factor functions O aij ___ >an (9i, . . . , 9 n ) by eq. (U) on page [| for the specific 
order of the rapidities 9\ > ■ ■ ■ > 9 n . For all other arrangements of the rapidities the 
functions Oq(9) are given by analytic continuation. Note that in general this analytic 
continuation does not provide the physical values of the form factor. These are given for 
ordered rapidities as indicated above and for other orders by the statistics of the particles, 
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of course. The Oa(9) are considered as the components of a co-vector valued function 
O\... n {0) £ V\... n — (X^ 1 "' 71 ) which may be depicted as 



01..M) = C ° ) • 

Now we formulate the main properties of form factors in terms of the functions C?i... n (0) 
which follow from general LSZ-assumptions and "maximal analyticity" . The later condi- 
tion means that Oi„, n (9) is a meromorphic function with respect to all 0's and all poles 
in the 'physical' strips < lm9ij < n (9^ = 9i — 9j i < j) are of physical origin, as for 
example bound state poles. 



Properties: The co-vector valued function O\... n {0) is meromorphic with respect to all 
variables 9\, . . . , 9 n and 

(i) it satisfies the symmetry property under the permutation of both, the variables 0j, 9j 
and the spaces i,j at the same time 

■ ■ ,di,9j, . . .) = 0...ji...(- ■ ■ , 9j, 9i, . . . ) Sij(9i — 9j) 

for all possible arrangements of the 

(ii) it satisfies the periodicity property under the cyclic permutation of the rapidity 
variables and spaces 

Ox... n (9i, 9 2 , ■ ■ ■ , 9 n , ) = 2 ...ni(9 2 , ■ ■ ■ , 9 n , 9i — 2m)a i 



{Hi) and it has poles determined by one-particle states in each sub-channel. In particular 
the function Oa(9_) has a pole at (9 12 = Z7r such that 

Res 0i.„„(0i, ...,9 n ) = 2i C 12 3 ... n (9 3 , . . . , 9 n ) (1 - S 2n . . . S 23 ) 

where C12 is the charge conjugation matrix. 

(iv) If the model also possesses bound states, the function Oa(9) has additional poles. 
If for instance the particles 1 and 2 form a bound state (12), there is a pole at 

(12) (12) 

9 12 = iu\ 2 (0 < u\2 < tt) such that 



Res^ 0i a ... n (0i, 9 2 , ... , n ) = ( i2)...„(0(i2), • ■ • , 9 n ) V2T^ 



(12) 

where the bound state intertwiner T 12 and the relations of the rapidities 61, 9 2) 9n 2 \ 



and the fusion angle u^f 1 [JTT 
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(v) Naturally, since we are dealing with relativistic quantum field theories we finally 
have Lorentz covariance 

0i... n (0i + fi, . . . , e n + n) = e s ^ c?i... n (0i, . . . , e n ) 

if the local operator transforms as O — > e s/1 where s is the "spin" of 0. 

In the formulae (i) the statistics of the particles is taken into account by S which 
means that — —Sn if both particles are fermions and S12 = S12 otherwise. In (ii) 
the statistics of the operator O is taken into account by oq\ = — 1 if both the operator 
O and particle 1 are fermionic and a 01 = 1 otherwise. 

The property (i) — (iv) may be depicted as 




We will now provide a constructive and systematic way of how to solve the properties 
(i) — (v) for the co-vector valued function / once the scattering matrix is given. These 
solutions are candidates of form factors. To capture the vectorial structure of the form 
factors we will employ the techniques of the algebraic Bethe Ansatz which we now explain 
briefly. 



2.3 The 'off-shell' Bethe Ansatz co-vectors 

As usual in the context of algebraic Bethe Ansatz we define the monodromy matrix as 
Ti... n fl(0, 9q) = S 10 (9 1 — 9 ) S 2 o(9 2 — 0o) - " S n0 (9 n — 9 ) = 



It is a matrix acting in the tensor product of the "quantum space" y x --- n = V± <g> • • ■ <g>V n 
and the "auxiliary space" Vo (all V{ = C 2 = soliton-anti-soliton space). The sub-matrices 
A, B, C, D with respect to the auxiliary space are defined by 











1 


2 


n 






71... n>0 (M 



Ai... n (£,z) Bi... n (£,z) 

C!...n(M £>!...„(£,*) 
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A Bethe Ansatz co- vector in V\... n is given by 

*l...n(£,*) = tt 1 ... n C 1 ... n (6,Z 1 )---C 1 ... n (6,Z r , 



0, 



0, 



9, 



(5) 



. , z m ). Usually one has the restriction 2m < n and the charge of the 
2m = number of solitons minus number of anti-solitons. The solitons 



where z_—\Z\ 
state is q = n 

are depicted by j or <— and anti-solitons by J. or — >. The co-vector fii... n is the "pseudo- 
vacuum" consisting only of solitons (highest weight states) 



tti...n =T ®---<8> T • 



It satisfies 



^l...n Bi_ n (9, z) 
^l...n^l...n(£, z) 

fil...»-Dl.„n(£, 





n 

n o(0< - ^)fii...„ 



n^-^)^i...n. 



The eigenvalues of the matrices A and .D, i.e. the functions a = —a and b = —b are given 
by the amplitudes of the scattering matrix @. In the following we use the co- vector 
^ i...n(9, z) in its 'off-shell' version which means that we do not fix the parameters z by 
means of Bethe Ansatz equations but we integrate over the z's. 



3 The general form factor formula 

In this section we present our main result. We derive a general formula in terms of an 
integral representation which allows to construct form factors i.e. matrix elements of local 
fields as given by eq. ([!]). More precisely, we construct co- vector valued functions which 
satisfy the properties (i) — (v) on page ||. 

As a candidate of a generalized form factor of a local operator O(0) we make the 
following Ansatz for the co- vector valued function 

£>i...n(0)= / dz x --- [ dz m h(6,z)p°(6,z)y 1 ... n (6,z) (6) 

JC 9 _ JC e _ 

with the Bethe Ansatz state ^i... n (@,z) defined by eq. (|5]). For all integration variables Zj 
(j = 1, . . . , m) the integration contours Cg consists of several pieces (see figure |l|): 

a) A line from — oo to oo avoiding all poles such that Im#j — 7r — e < Im Zj < Im 9i — n. 

b) Clock wise oriented circles around the poles (of the 0(0, — Zj)) at Zj — 9i (i — 
l,...,n). 
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o 6 x + m{2v-l) 

Q e n + m{2v-l) 6 2 + iix(2v-l) 

• Q\ + m 

• O n + in •9 2 + m 

fi + i-n[y - 1) 
o O n + w(y - 1) o 2 + Z7r(f - 1) n?\ e x 

(T)0 n ... $02 



oQn-i-K 



o 1 -in 



n — mv 2 — i'KV 



J 9i - mv 

• X - 2m 



• 0„ — 2iri • 2 — 2m 

Figure 1: The integration contour Cg_ (for the repulsive case v > 1). The bullets belong to 
poles of the integrand resulting from u(0i — Uj) 4>{0i — Uj) and the small open circles belong 
to poles originating from t(0i — Uj) and r{9{ — Uj). 

Let the scalar function (c.f. 

n m 

h(0,z)= n nynn^-^) n m 

l<i<j<n i=l j=l l<*<i<"i 



be given by 



and 



T{Z) = <P(Z)<P(- Z ) ' M = F(z) F(z + in) (8) 



F{0) = sm±0f£\0) 

ft\8) = c::p ^ {l ~ v)t l-«Mht(l-0/(i7r» 

Jo t sinh \vt cosh \t 2sinht 

The function F(0) is the soliton-soliton form factor fulfilling Watson's equations 

F(0)=F(-0)a(0) = F(2m-0) (10) 

with a(0) = —a(0) where a(0) is the soliton-soliton scattering amplitude. 

Remarks: 

• Using Watson's equations (|T0"| ) for F(z), crossing (|) and unitarity ([|) for the sine- 
Gordon amplitudes one derives the following identities for the scalar functions <j)(z) 
and t(z) from the definitions (^) 

1 a(z — 2-ri) 

4>{z) = (f)(tn -z) = j— <j>[z - m) = t— <p(z - 2m) , 

b(z) b{z) 
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t{z) = t(-z) = -±± -A— ( t(z - 2m) 

a(z) o(2irt — z) 

where b(z) is the soliton-anti-soliton scattering amplitude related to a(z) by crossing 
b(z) = a(m — z). 

The functions <fi(z) and t(z) are of the form 

const f 00 dt sinh |(1 - v)t (cosht(| - z/(m)) - l) 



— — exp / - =- — 

sinh z Jot smh | vt smh t 

t(z) = const, sinh z sinh zj v 

• The function h(6_,z) and the state ^i... n (^, z) are completely determined by the 
S-matrix. 

In contrast to the functions F(z), (j)(z) and t(z) the 'p-function' p (9, z) in the integral 
representation (|6]) depends on the local operator 0(x),in particular on the spin, the charge 
and the statistics. The number of the particles n and the number of integrations m are 
related by q = n — 2m where q is the charge of the operator 0{x). The p-functions is 
an entire function in the Zj (j = 1 , . . . , m) and in order that the form factor satisfies the 
properties (i) — (v) it has to satisfy the following 

Conditions: The p-function (6_, z) (where n is the number of particles and the number 
of variables 9) satisfies 

(*') Pnid'Z-) is symmetric with respect to the #'s and the z's. 

(W) Pn(9, z) = <JoiPn(- ■ ■ j @i ~ . . . , z) and it is a polynomial in e ±Zj (j = 1, . . . , m). 

The statistics factor o 0i is —1 if the operator 0(x) and the particle i are both 
fermionic and +1 otherwise. 

f Pn(0i = 9 n + in,~9,9 n ,~z,z m = 9 n ) = -p%_ 2 (9,z) +p^{9) 
(Hi') < _ m x ~ 

I P%(9i = 9 n + i7r,l9 n ;~z,z m = 9,) = cr m - P°- 2 (il) +P {2) (9) 
\ m 

where 9 = (9 2 , ■ ■ ■ , O n -i)i I = • • • z m-\) and where p^' 2 ^(£) are independent 
of the z's and non- vanishing only for charge less operators 0(x). The constant x 
depends on the coupling and is given by (see formula (B.7) in [|l|) 

/i 0) (o) N2 

x = — 



47T 



Uv' 



the bound state p-functions are investigated in section || 
v ') + + A 4 ) = eSti Pn(di£) where s is the 'spin' of the operator 0(x) 
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As an extension of theorem 4.1 in we prove the following theorem which allow to 
construct generalized form factors. 



Theorem 1 The co-vector valued function Oi,,, n (9) defined by jj^) fulfills the properties 
(i), (ii) and (in) on page if the functions F (9), <f>(z) andr(z) are given by definition (j^) 
- (j^) and if the p-function p® (9, z) satisfies the conditions (i') — (Hi'). 



This theorem is proven in [11 



Remarks: 



• The number of C-operators m depends on the charge q = n — 2m of the operator 
O, e.g. m — (n — l)/2 for the soliton field ip(x) with charge q — 1 and m = n/2 for 
charge- less operators like tpip or the energy momentum tensor T^ v . 

• Note that other sine-Gordon form factors can be calculated from the general formula 
(H) using the bound state formula (iv). 

• The general representation of form factors by formula (|6]) is not specific to the sine- 
Gordon model. It may be applied to all integrable quantum field theoretic model. 
The main difficulty is to solve the corresponding Bethe Ansatz. 



4 Examples of "p-functions" 

In this section we propose the p-functions for various local operators. Since the charge 
of the operators which we consider is zero the number of integrations m and the number 
of particles n are related by m = n/2 and form factors are non- vanishing only for even 
number of particles n — 2,4, ... . We consider p-functions of the form 

tm m \ 

P j(p")jy* + P °(p»)J2 e ~ LZi ( u ) 
j=i i=i / 

where P M is the total energy momentum vector of all particles. The integrals in (|6]) 
converge for L < (1/V + l)(n/2 — m + 1) + 1/u. For large values of L the form factors 
are defined in general as the analytic continuations of the integral representation from 



sufficiently small values of v to other values. Obviously the p-functions (p]) satisfy the 



conditions (%') — (Hi 1 ) on page ECU]. From the property (Hi') we obtain the recursion relation 



for the normalization constants 



, ^ r /n 1 

*n-2~ 



N ° = N °" ^ N? = N?-^*r- 1 



m mi 

1 'x) 



where iV^ follows from the two-particle form factors. For the local operators Af ^ 
J\f [ipj 5 ^ (x), the current j fl (x) = J\f [^7 M ^] (x), the energy momentum tensor T^ u (x 



11 



ipsyt 1 Q v ip ( x ) — g^ v C MT and the infinitely many higher conserved currents J^[x) we 



V 



propose the following p-functions 

tn m 
E e ^E e2l -E ee! E e ~" 
i=l i=l i=l 8=1 / 

/ n m n m 

Nt^ E e ~ ei E e "+E eei E e ~" 



pj 



V 



T+ - 



p 



p 



n m 



0,z) 
0,z) 

e,z) 



J=l i=l 
-1 



=1 i=l 

m n m 



E^'E^'+E^'E 1 

\i=l / \i=l i=l i=l i=l 

n / n \ ~ / n m n m 

< E « tt ( E eT "' ) E^'E «" -E e<>, E e ~ 



i=l 



t=l 



,i=l i=l 
n m 



1=1 i=l 



-jv? E^'E^'-E^'E^ 

\j=l i=l i=l i=l 

n m 

±N^J2 e±6i 12 eLZi > (i = ±l,±3,...). 



i=l i=l 



The identification with the operators has been made by comparing the exact results with 
Feynman graph expansions. Properties as charge, behavior under Lorentz transforma- 
tions etc. also become obvious [|llj] The fundamental sine-Gordon bose field tp(x) which 
correspond to the lowest breather is related to the current by Coleman's formula 

e^d uV = -jf or d ± V = ±jj ± . 

This implies for the p-function 

27T2 



n n 



-1 



n m 



n m 



pS(».2) = Jv^(E« , E e "') lE^E^+E^E 6 " 



i=i i=i 



,i=l i=l i=l i=l 



5 Soliton Breather form factors 



We calculate breather form factors starting with the general formula (^) for the soliton 
form factors using the property (vi) on page [| (for details see JTI], |I3|). The 6 fc -breather- 
(n — 2)-soliton form factor is obtained from Oi23...n(£) by means of the fusion procedure 



(iv) with the fusion angle given by u 12 
£ = 0(i 2) = i(0i + 2 ) and£' = 3 ,...,0, 



(12) 



li 



7r(l — fcv) the bound state rapidity 



Res O 



i 12 =iu 



123. ..nl 



0(l 2 )3...n(£,£)V^lf 2 2 



where the bound state intertwiner 11 . 
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For 9 12 — > ivP*' there will be pinchings of the integration contours in formula (|6|) at 
the poles Zi = z"' =92 — iirlv = £ — |«7r(l — fcf + 2lv) for Z = 0, . . . , k and i = 1, . . . , m. 
Using the pinching rule of contour integrals and the symmetry with respect to the m z- 
integrations we obtain 



k 

Res Oi23...n(9) = Res (— 2m) Res / dz® • • • / dz. 

_2=iu( k ) 12 =i u ( k ) zi=zW J Cg J Cg 



After a lengthy calculation [[II] we obtain for the case of the lowest breather the one- 
breather- (n — 2)-soliton form factor 

l 

3 ... n (Z,9!) = l[F sb (Z-9 l ) H F(9 ij ) Y^i-iyHpit-BiJ) 

2<i 2<i<j 1=0 2<i 



X 



/ dz 2 --- dz m JJxCC-^O 

J Cq_ JCe_ i^j 

xjinw-^) n <zii)p°^~,z {i \d)^..M,d) 



with z_' = (z2, . . . , z m ). The soliton-breather form factor has been introduced as 
F sb {9) = K sb {9) sin ±0 /i 0) (9) 

(0) r g cosh|z/t 1 - cosht(l - 9/ (in)) 

1st i ) ex py o t cosh | t 2sinht 

-cos 5(1- !/)/£(£(! -*)) 



K sb (9) 



sinh \{9 - f (1 + i/)) sinh i(0 + f (1 + i/)) 



The normalization has been chosen such that F sb (oo) = 1. The function -E(z/) was used 
in |5], m Also we have introduced the short notations 

sinh 2^ 



sinhi(e + f(l + (-l)V)) 



sinh±(£-f(l + (-l)'z/)) 
The following identities have been used 

F(9 1 -9 i )F(9 2 -9S(0i-z {l) ) = F sb (£ - 9 i )p(£ s - 9 h I) 



for 9\/2 — £ ± y(l — v), = £ — y(l — (—1)'^). The new p-function is obtained from 
the old one by 

f{S&,z®,d) = md{u)p°{i + \iu { - l \i- |iit (1) ,0',z (,) ,z') 
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where the constant d(v) is given by 



:-'\ sin \kv 



Iterating the procedure above we obtain the r-breather-s-soliton form factor with 
2r+s = n, the breather rapidities £ = . . . , £ r ) and the soliton rapidities # = (#1, . . . , S ) 



l<i<j<r i=l j=l l<«<i<s 

1 1 

X 



Z 1= / r =0 l<i<j<r V S ^ 7 



r s „ „ r m 

n n ^ ~ ^ z ) / dzr+i I dz m\\w x(ti - ^, o 



r s 

X 

i=l j=l " " i=l j=r+l 

s m 

i=l j=r+l r<i<j<m 

again with zf*' 1 = ^ — y(l — (— 1)^//), (i = 1, . . . , r). The two-breather form factor has 
been introduced as 

l t r(0), 



F bb (0 = K bb (0 sin ±U b Y(0 



(0)/m _ _ / ut o ws ^2 



oc 



dt o cosh(| - v )tl- cosht(l - 9 /(in)) 



fbb {d) ~ eXP 7 t " cosh |t 2 sinh/ 

/ — cos \-kvIE(v) 

K bb (9) = 



sinh \ (9 — mv) sinh |(0 + inv) 
The normalization has been chosen such that F bb (oo) = 1. It has been used that 

F sb (^ - 9 3 )F sb (^ - 9 4 ) - 9 3 , h)p{ii ~ 04, h)x{ti - 4 2 \ h) 



ti \ / , /, ism iris 

= F bb £ 12 )(l + (l 1 -l 2 — -— 
V smh£ 12 

for ^3/4 = £ 2 ± y(l — v). The new p-function is obtained from the old one by 

p&r,^,*") = ^H^p^ + i^),^-^ 1 ),...,^^ ,---,/) • 

In particular for n = 2r = 2m we get the pure lowest breather form factor 

o(e = n^««)E-"E(- 1 ) ,1+ -"- Hp n ( 12 ) 

i<j i!=0 I r =0 l=i<j V Smn ^y/ 

and the pure breather p-function 

P(i,zU) = r\ d r (u)p + £i - . . . , 4^, . . . ) . 
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6 The quantum sine-Gordon field equation 



The classical sine-Gordon model is given by the wave equation 

Oip(t, x) + — sin P(p(t, x) = 0. 

and two particle sine-Gordon S-matrix for the scattering of fundamental bosons (lowest 
breathers) jlBI 

sinh 9 + i sin nu 

= ■ ua ~ 

smh 9 — i sin ixv 

where 9 is the rapidity difference defined by pip 2 = m 2 cosh 9 and v is related to the 
coupling constant by v = (3 2 / (8n — (3 2 ). 

From the S-matrix off-shell quantities as arbitrary matrix elements of local operators 
are obtained by means of the "form factor program" ||. In particular we provide exact 
expressions for all matrix elements of all powers of the fundamental bose field <p(t,x) 
and its exponential jVexp i^f(p(t, x) for arbitrary 7. Here and in the following M denotes 
normal ordering with respect to the physical vacuum which means in particular for the 
vacuum expectation value (0 | Af expij<p(t, x) \ 0) = 1. For the exceptional value 7 = j3 
we find that the operator □ _1 A/"sin (3ip(t, x) is local. Moreover the quantum sine-Gordon 
field equation^ 

n<p(t, x) + -AT sin p<p(t, x) = (13) 

is fulfilled for all matrix elements, if the "bare" mass y/a is related to the renormalized 
mass by0 

a = m 1 (14) 

sin ixv 

where m is the physical mass of the fundamental boson. The factor . wv modifies the 

1 J sin -kv 

classical equation and has to be considered as a quantum correction. For the sinh-Gordon 
model an analogous quantum field equation has been obtained in f29|Q. Note that in 
particular at the 'free fermion point' v — > 1 (/? 2 — > An) this factor diverges, a phenomenon 



which is to be expected by short distance investigations |18[ . For fixed bare mass square 
a and v — > 2,3,4, .. . the physical mass goes to zero. These values of the coupling are 
known to be specific: 1. the Bethe Ansatz vacuum in the language of the massive Thirring 



model shows phase transitions [HJ and 2. the model at these points is related [f2(], |21], ^2 
to Baxters RSOS-models which correspond to minimal conformal models with central 
charge c = 1 — 6/ (u(u + 1)) (see also pgji ). 



3 In the framework of constructive quantum field theory quantum field equations where considered in 

00- __ 

4 Before such a formula was found in |25|, |26| . 

5 It should be obtained from (|lh by the replacement (3 — > ig. However the relation between the bare 
and the renormalized mass in j29fl differs from the analytic continuation of ( |l4| ) by a factor which is 
l + 0(/3 4 )^f. 
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Also we calculate all matrix elements of all higher local currents J^(t,x) (M = 
±1, ±3, . . . ) fulfilling d^J^itjx) = which is characteristic for integrable models. The 
higher charges fulfill the eigenvalue equation 



dxJ° M (x) - {pf) M ) I Pi, ■ ■ ■ ,PnY n = 0. (15) 

i=l 



In particular for M = ±1 the currents yield the energy momentum tensor T^ v = T Vii with 



d^T^ = 0. We find that its trace fulfills 



T^(t,x) = -2^ (l - §- ) (tfcoa /3<p(t,x) - 1) . (16) 



(3 Z V 8^ 

This formula is consistent with renormalization group arguments [ H| , |23| . In particular 
this means that f3 2 /An is the anomalous dimension of cos f3(p. Again this operator equation 
is modified by a quantum correction (1 — (3 2 /8tt). Obviously for fixed bare mass square 
a and (3 2 — > 8tt the model will be conformal invariant which is related to a Berezinski- 
Kosterlitz-Thouless phase transition | 24j . The proofs of the statements (|13D - (p!6|) is 
sketched in the following together with some checks in perturbation theory. The complete 



proofs will be published elsewhere [J I 



A form factor of n fundamental bosons (lowest breathers) is of the form || 

fSiB = N°K°{6) n F{0ii) 

l<i<j<n 

where N® is a normalization constant, dij = 6i — 6j and F{9) is the two particle form 
factor function. It fulfills Watson's equations 

F(6) = F(-6)S(6) = F(2m - 6) 

with the S-matrix given above. Explicitly it is given by the integral representation 

m = N r dt (cosh \t - cosh(l + u)t) (1 - cosht(l - p) 
J Q t cosh |t sinh t 

normalized such that F(oo) = 1. The K-function K®{6) is meromorphic, symmetric and 
periodic (under 9i — > 6i + 2m). 
From (O) we have 



nfffl-E-EH)"^ n (i+c-ySf) ! '?(«.d (17) 

l 1= l n =0 l<i<j<n ^ % 'j ' 



where Z{ = 9i — y (1 + (2k — l)^). The dependence on the operator is encoded in the 
'p-function' p®. It is separately symmetric with respect to the variables 9 and z and has 
to fulfill some simple conditions in order that the form factor function f® fulfill some 
properties |j| ||. These properties follow (see 0) from general LSZ-assumptions and in 
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additions specific features typical for integrable field theories. In particular the recursion 
relation holds 

Res ff{9 Xt . . . , 6 n ) = 2i /„°_ 2 (fl 3 , ■ ■ ■ , 6 n ) (1 - S(6 2n ) . . . S(9 23 )) . (18) 

012=i7T 

Here we will not provide more details but only give some examples of operators and their 
corresponding p-functions: 

1. The correspondence of exponentials of the field and their p-function[] is 

n 

WV W «-> Y[ e {2h - 1)iwu ^ /l3 (19) 
i=i 

for an arbitrary constant 7. 

2. Taking derivatives of this formula with respect to 7 we get for the field and its 
powers 

Ny N <-> fe( 2 ^ - X ) 1 • ( 2 °) 

3. Higher currents (for M = ±1, ±3, . . . ) correspond to the p-functions 



for n = even and zero for n = odd. For M = ±1 we get the light cone components 
of the energy momentum tensor T pcr = J? with p, a — ± (see also [p9|). 

In order to prove equations as for example (|13|) and (|T^) we consider the corresponding 
p-functions and their K-functions defined by ([T7j) . The K-functions are rational functions 
of the Xi = e d \ We analyze the poles and the asymptotic behavior and find identities 
by using induction and Liouville's theorem. Transforming these identities to the corre- 
sponding form factors one finds the field equation (|i~3"D and the trace equation fll6D up to 
normalizations. 

Normalization constants are obtained in the various cases by the following observa- 
tions: 

a) The normalization a field annihilating a one-particle state is given by the vacuum 
one-particle matrix element, in particular for the fundamental bose field one has 

(0\<p(0)\p) = \fz* 

6 For the sinh-Gordon model an analogous representation as ( |l7| ) together with this p-function was 
obtained in |3Q] by different methods. 
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with the finite wave function renormalization constant calculated in || as 

-v ( i r v t 

Zf = (l + u)^^esxpl— I —dt 



sin \ 7T J smt 

where Z^ = 1 + 0(/3 4 ). This gives the normalization constant 

n[ 1] = Vz^/2 (21) 

for the form factors of the fundamental bose field and which are obtained from the 
p-function of (f20|) for N = 1. 



b) If a local operator is connected to an observable like a charge Q = J dxO(x) we use 
the relation 

(p'\Q\p) = q(p'\ p)- 

For example for the higher charges we obtain 

iM m M+i j 02 gin7rz/ 



*i u = n . m - - with ^- = Z- 



2 sinnvsm If n vF (in) F(m) 



inv iris 



c) We use Weinberg's power counting theorem for bosonic Feynman graphs pTTflf 7 ] . For 
the exponentials of the boson field O = Me %1Lp this yields in particular the asymp- 
totic behavior 

fiiOifi*, ...) = f?{e 1 ) f°_ x {e 2 , ...) + o(e~^) 

as Re 6*i — > oo in any order of perturbation theory. This behavior is also assumed 
to hold for the exact form factors. Applying this formula iteratively we obtain from 
( |PT|) with ( |l9|) the following relation for the normalization constants of the operators 

N2 = (Nl) n (n>l). 

d) The recursion relation (pl|) relates N n+2 and N n . For all p-functions mentioned above 
we obtain 

N n+2 = N n - 2 — — (n>l). 

Sin TTVr {17T) 

Using c) and d) we calculate the normalization constants for the exponential of the 
field A/V w and obtain 

iV7 = VZ^—. (22) 
1 2txv K J 



The normalization constants (|2lD and ( \22\) now yield the field equation (|13| ) with the 
mass relation (Jl4j) . The statement ( fl6|) is proved similarly. The eigen value equation (JTSj) 



18 




(a) (b) 
Figure 2: Feynman graphs 

is obtained by taking the scalar products with arbitrary out-states and by using a general 



crossing formula [11 



All the results may be checked in perturbation theory by Feynman graph expansions. 
In particular in lowest order the relation between the bare and the renormalized mass 
([14]) is given by Figure 1 (a). It had already been calculated in || and yields 



m * = a [l-l(£)\o^ 

which agrees with the exact formula above. Similarly we check the quantum corrections 
of the trace of the energy momentum tensor ( |T6| ) by calculating the Feynman graph of 
Figure 1 (b) with the result again taken from || as 

(p\Mcoap<p-l\p) = -P 2 (l + Q + 0(/? 6 ) . 

This again agrees with the exact formula above since the normalization given by eq. (p~5|) 
implies {p\T^\p) = 2m 2 . All other equations have also been checked in perturbation 



theory flT3|. 
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